The electronic structure of Nalr03, Mott insulator or band insulator? 



Liang Du, Xianlei Sheng, Hongming Weng0 and Xi Da{t] 

Beijing National Laboratory for Condensed Matter Physics, 
and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 

(Dated: October 11, 2012) 

Motivated by the unveiled complexity of nonmagnetic insulating behavior in pentavalent post- 
perovskite NalrOs, we have studied its electronic structure and phase diagram in the plane of 
Coulomb repulsive interaction and spin-orbit coupling (SOC) by using the newly developed local 
density approximation plus Gutzwiller method. Our theoretical study proposes the metal- insulator 
transition can be generated by two different physical pictures: renormalized band insulator or Mott 
insulator regime. For the realistic material parameters in NalrOs, Coulomb interaction U — 2.0(J = 
t//4) eV and SOC strength 77 = 0.33 eV, it tends to favor the renormalized band insulator picture 
as revealed by our study. 

PACS numbers: 71.15.Mb, 71.27.+a, 71.30.+h 



Recently there are increasing research activities on the 
spin-orbit coupling (SOC) driven metal-insulator transi- 
tion (MIT) in 5d transition metal compounds.^ ^ Com- 
pared with 3d and 4d orbitals, the 5d ones are spatially 
more extended, leading to smaller on-site Coulomb re- 
pulsive interaction. However, since the effect of Coulomb 
interaction is comparable with SOC, the interplay be- 
tween these two factors can lead to unexpected inter- 
esting phenomena in these materials.^ The first well 
studied material of this type is Sr2lr04,^ where t2g and 
Cg orbitals are separated by large crystal-field with the 
lower t2g orbitals filled by five electrons. Due to strong 
SOC in the system, the t2g orbitals are split into nearly 
full- filled fourfold jeff = 3/2 states and half- filled twofold 
Jeff = 1/2 states. As a result, the system can be simplified 
to an effective one-band half-filled system with reduced 
bandwidth, which has Mott insulator ground state with 
antiferromagnetic long range order. Further theoretical 
studies show that both SOC and on-site Coulomb interac- 
tion are essential to explain the Mott insulator behavior 
in this material, which has also been supported by recent 
experimental studies including angle-resolved photoemis- 
sion spectroscopy, optical conductivity, x-ray absorption^ 
and resonant x-ray scatting^ measurements. 

Recently, another novel post-perovskite compound 
NalrOs with pentavalent iridium Ir^+ ions was synthe- 
sized by R. Cava's group in Princeton. The transport 
and susceptibility measurements show that NalrOs is an 
insulator without magnetic order. However, electronic 
structure calculations based on density functional theory 
(DFT) with local density approximation (LDA) found 
significant density of states (DOS) at the Fermi level, in- 
dicating the material to be metallic, which is in strong 
contrast with the experimental observations. The fail- 
ure of DFT type calculation on predicting the basic elec- 
tronic structure of the material implies that the strong 
correlation effect may play an important role here, which 
is poorly treated by LDA alone. The DFT calculation 
shows that in NalrOs the energy bands near the Fermi 
level are mainly formed by the Ir t2g orbitals from its bd 
shell. The SOC splits the t2g bands into two groups of 



bands with effective total angular momentum jeff = 1/2 
and Jeff = 3/2 respectively. Since for NalrOs there are 
totally four electrons remaining in the t2g bands, large 
enough SOC will naturally lead to band insulator phase 
with fully occupied jeff = 3/2 bands and empty jeff = 1/2 
bands. While the LDA calculation predicts a metallic 
phase, simply because the spin-orbit splitting is still sev- 
eral times smaller than that of the t2g bandwidth and is 
not strong enough to generate a band insulator by SOC 
alone. 

In this type of system, there are two possible ways 
for the correlation effect induced by local Coulomb inter- 
action to generate insulating behavior. The first one is 
called "renormalized band insulator", in which the cor- 
relation effect reduces the effective bandwidth and en- 
hances the effect of SOC leading to a band insulator 
phase with "renormalized" band structure. The second 
possibility is the Mott insulator with completely van- 
ishing of quasiparticles, which is caused by strong local 
Coulomb repulsive interaction. Unlike the Mott insula- 
tor phase in Sr2lr04, which has five bd electrons on each 
Ir ion leading to magnetically ordered phase at zero tem- 
perature, in NalrOs there are four electrons in the t2g or- 
bitals, which leads to nonmagnetic atomic ground states 
with spin and orbital moments canceling each other even 
in the strong coupling limit with very large Coulomb in- 
teraction and weak SOC. Since there is no symmetry dif- 
ference between the above two possible insulator phases, 
these two "phases" should be adiabatically connected to 
each other and there will be no phase transition but 
crossover between them in the parameter space. While 
it is then interesting to ask that for this particular ma- 
terial NalrOs, is it more close to a "renormalized band 
insulator" or "Mott insulator"? 

In the present letter, we apply the newly developed 
LDA+Gutzwillei^^ method to study the interplay be- 
tween SOC and local Coulomb interaction among bd elec- 
trons in NalrOs- We find that both the above two effects 
are important to explain the insulating behavior in this 
material and this particular material would be better de- 
scribed by the picture of "renormalized band insulators" 
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than that of Mott insulator. 




FIG. 1: (Color online) The experimental crystal structure of 
NalrOs is shown in the left panel. The green (large) ball 
represents Na atom. Red (small) and darkslategray (middle) 
balls represent O and Ir atoms, respectively. Total density of 
states and the projected one for Ir d orbitals obtained from 
LDA calculation of NalrOs are shown in the right panel. 

Firstly, we carry out the electronic structure calcula- 
tion by DFT with LDA using the experimentally deter- 
mined post-perovskite crystal structure/^ which is shown 
in Fig{l] In such structure each Ir atom is surrounded 
by six oxygen to form an octahedron. These octahedrons 
are connected by sharing Conner oxygen along c-axis and 
by sharing edge oxygens along a-axis, which forms IrOa 
sheet stacking along 6- axis separated by Na ion layer. 
Our LDA calculation indicates that the coupling between 
neighboring IrOs sheets are quite weak, being consistent 
with the two-dimensional electron transporting property 
found in transport experiments. To search for the pos- 
sible magnetic ordering of Ir ions, four types of mag- 
netic orderings, namely, ferromagnetic ordering in both 
a- and c-axis (noted as aFcF), ferromagnetic ordering 
along a- while anti-ferromagnetic ordering along c-axis 
(noted as aFcA), anti- ferromagnetic ordering along a- 
and ferromagnetic along c-axis (noted as aAcF) , and anti- 
ferromagnetic ordering along both a- and c-axis (noted 
as aAcA), are studied. The first-principles calculations 
of these states are performed by using full potential all 
electron method implemented in WIEN2k software pack- 
age in order to get highly accurate total energies. The 
LDA calculations can get only nonmagnetic solution and 
density of states in Figfllshow that Ir 5d t2g and Cg or- 
bitals are well separated by about 3.5 eV under octahe- 
dral crystal-field. 

The bandwidth of t2g orbitals are about 3.0 eV, which 
is quite large since the direct overlap of neighboring t2g 
orbitals is possible due to the edge-sharing connection 
along a-axis. The main physics is dominated by Ir t2g 
orbitals around the Fermi level. Since both the correla- 
tion effect and the SOC are important in 5d transition 
metal compounds, we further apply the LDA+U+SOC 



to study the electronic structure of this material with 
U varying from 1.0 to 7.0 eV. Our numerical results 
show that the LDA+U+SOC calculations always con- 
verge to nonmagnetic and metallic solution even when U 
is as large as 7.0 eV. These results indicate the failure 
of mean field treatment of such a correlated system and 
profound infiuence of SOC on the electronic structures. 
In order to well describe the observed insulating behavior 
in NalrOa, we have further performed LDA+Gutzwiller 
calculations combining the DFT with the Gutzwiller vari- 
ational method, which can treat the correlation effects 
more precisely. 

We firstly construct an accurate low energy model 
Hamiltonian, which can well catch the crystal-field split- 
ting, the hopping parameters among active t2g local or- 
bitals and the SOC in real material. To do this, the 
projected atomic Wannier functions^^ (PAW) are con- 
structed for the t2a orbitals of Ir ion by using OpenMX 
software packag#^ within LDA calculation. This ap- 
proach has been used to treat 3d transition metal t2g 
orbitals successfully.'^ In one unit cell, there are two 
non-equivalent Ir ions and the t2g orbitals are defined 
in each local coordinates. The obtained PAW orbitals 
are quite localized and the band structure obtained by 
first-principles calculations can be well reproduced by the 
tight-binding Hamiltonian using these basis. It is shown 
that the dxy-\ike PAW is lower than the nearly degen- 
erated dxz- and dyz-\ike PAWs by about 0.52 eV, which 
is known as tetragonal crystal-field. The atomic SOC 
is added to the above tight-binding (TB) Hamiltonian 
and the effective SOC strength 77 is found to be 0.33 eV 
when fitted with the LDA+SOC calculations. As a re- 
sult, the effect of SOC splits jeff = 1/2 (higher in energy) 
and Jeff = 3/2 states by about 3r]/2 0.5 eV. There- 
fore, the strength of SOC and tetragonal crystal-field are 
comparable and they compete against each other in con- 
sidering of the orbital degeneracy. Implemented with the 
local Coulomb interaction terms among the t2g orbitals, 
the total Hamiltonian can be written as, 

H = Ht^Hu^H^^HA. (1) 

The first term describes the hopping process of electrons 
between local spin-orbitals "acr" and "6cr'", 

where a denotes electronic spin, and a represents the 
three t2g orbitals with a = 1,2,3 corresponding to 
dyz^dzxidxy orbitals respectively. The rest terms of the 
Hamiltonian are all local terms expressed by Hl^^ = 
Hl^ -\- Hl^-\- i^^, which contains Coulomb interaction i^* , 
SOC JY* and tetragonal crystal-field splitting i^^. (In 
the following, the site index is suppressed for sake of sim- 
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plicity). 

a a<b,aa' a<b,a 

(4t^«i4i^^t + dl^dl^db^dbi + h.c}j , (3) 

a<b 

Hrj = r]{aa\l + ks^\ba')di^dba', (4) 

aa,ba' 

Ha= ^aaM'^ia^ba', (5) 

aa^ba' 

where U (U') is the strength of intra-orbital (interor- 
bital) Coulomb interaction and J describes the Hund's 
rule coupling. [/, U' and J satisfy Kanamori constraint 
U = -\-2J. I and s represent the orbital and spin angu- 
lar momentum operators, respectively. 77 represents the 
SOC strength and A describes the crystal-field splitting 
between the three t2g orbit als. We emphasize that in the 
present study we treat the Hund's rule coupling terms in 
a full rotational invariant way, which is crucial to guar- 
antee that the locking of spin and orbital spaces is purely 
due to SOC but not the Sz-Sz Hund's coupling term with 
artificially chosen z-direction. We also fix the ratio be- 
tween Hund's coupling J and Hubbard interaction U to 
be 1/4 throughout the entire paper. 

Next, we briefly introduce the Gutzwiller wave func- 
tion (GWF) used in this paper. The generalized GWF 
I^g) with rotational invariant local interaction terms can 
be constructed by acting a many-particle projection op- 
erator V on the uncorrelated wave function l^o), 

I^g) =7^|^o), (6) 

with 

P = n = n E ^ Wrr' Ir, R) r, R| . (7) 

R R rr' 

where |^o) is a normalized uncorrelated wave function 
in which Wick's theorem holds, |r,R) represents atomic 
eigenstates on site R and A(R)rr' are Gutzwiller varia- 
tional parameters to be determined by variational princi- 
ple. In our work, |r,R) are eigenstates of atomic Hamil- 
tonian Hioc- The expectational value of hopping terms 
Ht in our Hamiltonian can be expressed as: 

ij a (3 57 

(8) 

with (5, 7) being combined spin-orbital index and 



^ir = {I\v\r)^/W^\IWW^. (10) 



where |/) stand for the many-body Fock states and 
= (^0 1^7 1^0) with rij being occupation number op- 
erator for 7 state. The Gutzwiller variational param- 
eters A(R)rr' are determined by minimizing the above 
total energy. The quasiparticle weights in the general- 
ized Gutzwiller method is deflned as the eigenvalues of 
the Hermite matrix IZ^IZ. The detail numerical proce- 
dure for the rotational invariant Gutzwiller method can 
be found in reference [19H26]. The complete phase di- 
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FIG. 2: (Color online) The phase diagram in the plane of 
Coulomb interaction and spin-orbit coupling (SOC). There 
exists three different regions: Mott insulator, Band insulator 
and Metal. The blue star locating in band insulating region 
denotes the realistic parameters of NalrOs with Coulomb re- 
pulsive interaction JJ — 2.0 eV and SOC strength r\ — 0.33 
eV. The Hund's rule coupling J is fixed as [7/4. 
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FIG. 3: (Color online) The quasiparticle weight 2; as a func- 
tion Coulomb interaction for different SOC strength 77. Note 
here the orbital is selected to be the one with smallest quasi- 
particle weight. For r\ — 0.1(0.2) eV, there exists a first order 
transition to Mott insulator around U — 2.8(3.2) eV. While 
for r\ — 0.3(0.4) eV, the transition to band insulating state is 
around IJ — 2.2(1.2) eV. The Hund's rule coupling J is fixed 
as [7/4. 

agram in the parameter space spanned by SOC r] and 
Hubbard repulsive interaction U with fixed J jU = 1/4 
has been plotted in Figure|2] According to the strength 
of SOC, the whole phase diagram can be divided into two 
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FIG. 4: (Color online) At realistic SOC strength rj = 0.33 
eV, the band structure of NalrOs with [7=1.0 eV (left panel) 
and 2.0 eV (right panel) are calculated. For [/ = 2.0 eV case, 
clearly there exists a gap of about 500 meV. The Hund's rule 
coupling J is fixed as [7/4. 

regions. In the left part of the phase diagram, the SOC 
is much weaker than the tetragonal crystal-field, which 
splits the t2g orbitals into two-fold (including spin) low 
lying and four-fold high lying energy levels, respectively. 
With the lower two-fold bands being fully occupied, the 
rest of the system can be considered effectively as two- 
band system filled by two electrons. A typical Mott tran- 
sition happens when the Hubbard repulsive interaction U 
increases over the critical value, which is around 2.7 eV 
in the weak SOC limit. The quasiparticle weight 2; as a 
function of U for the weak SOC strength {r] = 0.1,0.2 
eV) is plotted in Figjsj which shows typical Mott transi- 
tions with strong first order nature. The behavior of the 
quasiparticle weight, both the concave shape of the curve 
and the first order nature for the Mott transition, agrees 
quite well with results obtained by studying the model 
Hamiltonians with semi-circle like density of states. The 
increment of SOC will further split the above two bands 
and transfer part of the electrons from one band to an- 
other, which will suppress the orbital correlation in the 
effective two-band system and raise the critical Uc for the 
system. The same increment of Uc as the function of en- 
ergy level splitting in the two-band Hubbard model has 
been also reported and discussed in detail by P. Werner 
et al.^2 

On the right part of the phase diagram, compared to 
the tetragonal crystal- field, the SOC becomes more dom- 



inant and we find a completely deferent behavior as we 
increase the interaction strength U. The quasiparticle 
weight z in the large SOC region (77 = 0.3, 0.4 eV) is 
again plotted as the function of U in Figj3| which in- 
dicates that the quasiparticles with non-zero weight z 
persists all the way up to above 4.0 eV. Detailed analy- 
sis of the corresponding band structure, which is plotted 
in Figj4] for SOC strength r] = 0.33 eV, indicates that 
there is a transition to band insulator around Uc = 1.8 
eV, above which a clear band gap appears between the 
forth and fifth bands as illustrated in Figj4j Unlike the 
Mott transition we discussed in the previous paragraph, 
the transition in the right part of our phase diagram is a 
typical Lifshitz transition characterized by the vanishing 
of the electron and hole fermi surfaces right at the tran- 
sition point. Therefore for this particular material, with 
the different strength of SOC, the on-site Coulomb in- 
teraction can induce Mott transition for weak SOC and 
Lifshitz transition for strong SOC cases. The realistic 
SOC strength of NalrOa has been fitted to be around 
0.33 eV, which is indicated in Figj2] by a blue star and 
clearly located at the right side of the phase diagram. 

In conclusion, our LDA+Gutzwiller calculation con- 
firms that NalrOs is a band insulator with renormalized 
bandwidth induced by the correlation effect among the 
5d orbitals. As shown in FigjSj the quasiparticle weight 
drops to about 0.65 ai U = 2.0 eV, which completely 
delimitates the overlap between the conduction and va- 
lence band leading to a band insulator with renormalized 
bandwidth. The SOC in this material plays a very impor- 
tant role, which leads to very different behavior for weak 
and strong SOC. In the former case, the on-site Coulomb 
interaction among the 5d orbitals will generate a Mott 
transition, after which the system is well described by lo- 
cal moments. While in the latter case, when the strength 
of SOC overwhelms the tetragonal crystal-field, the cor- 
relation effect will always generate band insulators with 
renormalized bandwidth. The correlation effect in these 
renormalized band insulator phase may also manifest it- 
self in the dynamical properties, i.e. the unusual exciton 
behavior, which will be discussed in our further publica- 
tions. 
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